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Based on the dielectric formalism in the generalized random phase approximation, we generalize 
the description of a Bose condensed gas to allow for a relative velocity between the superfluid and 
normal fluid. In this model, we determine the critical velocity dynamically as the transition point 
between stable and unstable dynamics. Unlike the zero temperature case, at finite temperature 
the relative critical velocity of a dilute Bose gas is lower than the sound velocity. This result 
illustrates one relevant difference that exists between a conserving and gapless approximation and 
other approaches. 
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One of the most fascinating aspects of the observed 
Bose condensation phenomenon is the superfluidity prop- 
erty i.e. the ability of the superfluid fraction of the Bose 
gas to move without apparent friction. No dissipation oc- 
curs provided that the relative motion of the superfluid 
with an obstacle, like an impurity or a local potential or a 
wall, does not exceed the sound velocity. Otherwise, the 
Landau criterion predicts that a phonon can be emitted 
decreasing the velocity of the superfluid pj . The situa- 
tion is drastically different if the obstacle corresponds to 
the normal component of the gas. In principle, the Lan- 
dau criterion cannot be applied to the normal fluid as 
a whole but to each of its individual thermal excitation. 
Unfortunately, some of them are highly energetic and 
display a relative velocity higher than the sound velocity. 
In order to explain the persistence of the relative motion 
between the normal and superfluid, a full microscopic 
description is necessary. A quantum kinetic equation ex- 
plaining this persistence has been proposed in |2| and has 
been subsequently extended in 0, 0, |j| . At low temper- 
ature i.e. ksT < gn (where g = Aira/m, a the scattering 
length and n the total gas density), this equation is an 
extension of the Beliaev approach for a non equilibrium 
Bose gas and predicts that the critical and sound veloc- 
ity are identical. However, the extension of the Beliaev 
theory including the presence of a thermal component is 
not conserving [f| . In particular the conservation law for 
the total mass density is missing in the systematic calcu- 
lation made in |5( using the Keldysh formalism. On the 
contrary, at high temperatures (ksT > gn), the conser- 
vation laws are fulfilled but at the price of having a gap 
in the energy spectrum rendering possible the friction- 
less motion 0i y, Of course, by the Hugenholtz-Pines 
theorem, a gap is forbidden and cannot be used to ex- 
plain the superfluidity phenomenon. As a consequence, 
a gapless and conserving approach is needed in order to 
improve the description in a large range of temperature. 
Recently, such an approach has been developed in both 
equilibrium and in non equilibrium 8] using the gen- 
eralized random phase approximation and has been suc- 
cessful to explain the low-lying excitation observed in a 
trapped Bose gas ||- 

In this letter, based on this more coherent approach, we 



delimitate the region for which a relative velocity between 
the normal fluid and the superfluid is possible at finite 
temperature. Unlike for other obstacles, the transition 
region from metastability to instability, due to counter- 
flow between normal and superfluid components, occurs 
for velocity lower than the sound velocity. The persis- 
tence is explained both by the equilibrium and non equi- 
librium approaches. In equilibrium, we introduce an 77- 
ensemble for which the condensed fraction evolves with 
a velocity v s . The study of the dynamical fluctuations 
around this equilibrium allows to establish a critical value 
of v s above which these become instable. In non equi- 
librium, a previous work has shown that the supression 
of the binary collision process between condensed and 
non condensed atoms allows precisely this persistence § . 
Above the critical velocity, this collisionless regime dis- 
appears leading to a damping of the relative velocity. We 
will show that the equilibrium formalism provides iden- 
tical results to the non equilibrium one, in particular for 
the dispersion relation of the collective excitations. 

To start with, we consider a homogeneous Bose con- 
densed gas populated with atoms with mass m in a vol- 
ume V. For any momentum k, n c is the density of the 
condensate and the density of the normal cloud in 
the mode k. We extend the formalism developed in 
to the case where the N c — Vn c condensate particles 
are moving with a momentum k s . For that purpose, we 
break the U(l) symmetry introducing the 77— ensemble by 
adding the term proportional to the field operator Ck s to 
the total many body Hamiltonian H. Thus, any thermal 
averaging is written as: 

(A) = T±(Ae-P [H - vc ^- v * c l> ] )/T±(e-P [H - nCk »- ri * c L» ] )(l) 

In this way, the condensate wave function (ck s ) = \^N C 
is adjusted to move with a velocity v s = k s /m. The 
field operator in the momentum space can be written as 
Ck = 5c k + (5 k ,k 5 \/A 7 c- 

The Green's functions descri bing the non condensed 
fraction can be written as usual jlfjj : 

G a0 (k,u) = - [ dre iUT (T T [Sc ka (T)ScU) (2) 
Jo 
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where <5c kct is <5c k for a — 1 and <5<4 k _ k for a = 2. 
The only difference with the previous approach jTJ is that 
the propagator line of momentum k is coupled to the 
one with momentum 2k s — k (see Fig.f). The Green's 
functions satisfy the Bcliacv-Dyson relations: 



G a p — G a p + G^Hy/jGa^ 



(3) 



where (k, uj) are the usual self-energy matrix and 
where the Green's functions of the non interacting system 
in presence of the momentum k s is given by: 



The main difference to is that the susceptibility func- 
tion is expressed in terms of the momentum q — k — k s 
relative to the condensate momentum while the Green's 
function is expressed in terms of the absolute momentum 
k. As a consequence, the momentum of the incoming and 
outgoing lines in the vertex function must be relabelled 
adequately as shown in Fig. (2). 

Once the regular functions are known, any physical 
quantities can be calculated. In the random phase ap- 
proximation, these functions are found to be: 
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where e k = k 2 /(2m). Solving we obtain: 
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where 

D = (w-q.k s /m-A_) 2 -(e q + e ks -At + ^+) 2 
+ S12S21 (6) 

A± = (£n ± £ 2 2)/2 and q = k k s . The susceptibility 
function is defined as: 



Ai r) (q,c) 



= Xnn(q»w)/(1 - SX°n(q>«>)) (12) 
= jN c /{l-g X l n M) (13) 

= g(2n - n c )S a ./3 

+ .9 2 ^ c xl(q,^)/(i-.9x°„(q^)) (14) 



These are expressed in terms of the susceptibility func- 
tion of the non interacting normal fluid: 



xL(q» w ) = z_ 
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where p q = Ek 4+q Ck - H can b e expressed in terms of 
the proper dielectric function \ n n through: 



Xm*(q, w) = Xrm(q, ^)/(l - 5Xrm(q, w)) 



(8) 



The proper part of a quantity is the sum of only those 
of its diagrams which remain connected after cutting a 
single interaction line. All the unknown functions can be 
expressed in terms of regular (i.e. proper and one-particle 
irreducible) quantities, denoted by an upper index: 

X nn (q,Lu) = x&2(q,u) + AM(q,w)G Q/3 (k,w)Aj J r) (q,w) 

(9) 
(10) 

(11) 
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Below the condensation point, the chemical potential 
is determined through the tadpole condition or the 
Hugenholtz-Pines theorem /i = e ks + <?(2n.— n c ) and con- 
tains a kinetic energy term of the condensate particle 
e ks . The particle momentum distribution for a normal 
fluid at rest reads in the Hartree-Fock approximation as 
K = l/[e^ k2 / 2m+29 "-^ - 1]. The equation of state is: 



7U/2(^- 



2gn) 



(16) 



where g n (x) = Y^jLi xj /j n an< 4 ^th = \j2-K/mk B T de- 
notes the thermal wavelength. The critical temperature 
is given by k B T c = 2nm(n/((3/2)) 2 / 3 . 

With all these building blocks modified in presence of 
k s , we get the results: 
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G u (k,w) 
Gia(k,w) 
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A(q,w) 



(a; - k s .q/m + e q )[l - 2gx°„(q,^)] +g»e[l + 2gx° TI (q,^)] 

A(q,w) 

gn c [l + 2g X ° nn (q,uj)} 
A(q,w) 

(w - k s .q/m) 2 - el X°„(q> w ) + 2n c e q [l + 5X°n(q, 



A(q,w) 

(w - k s .q/m) 2 - e 2 . [i _ 2 5 x°„(q, w)] - 2<?n c e q [l + 2srx°„(q ; "0 



(17) 
(18) 

(19) 
(20) 
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FIG. 1: Representation of the normal and anomalous Green's 
functions in presence of a condensate velocity k s 
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FIG. 2: Representation of the regular vertex part ' m pres- 
ence of a condensate velocity k s coupling a particle line of 
momentum q + k s with a condensate line with momentum k s 
and an interaction line with momentum q. 

We are looking for the zeroes of A(q, u>) the dielectric 
function. The real part and imaginary part of the roots 
to correspond to the collective excitation and the Landau 
damping respectively and have been determined previ- 
ously in the case of thermodynamic equilibrium 0, @- 
For low momentum excitations, we consider the limit 
u) — ► and |q| — > and introduce the complex veloc- 
ity of sound by w/|q| = c. The only difference is that 
now we must consider terms of the form k s .q/|q|. If the 
wavevector q is parallel to the critical velocity, we finnd 
the following generalization: 



(c-v s ) 2 + c| 



where 



1 f c 
Xn(s) = —j= \ 
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I e t 2 -0(v-2gn) _ I 



(21) 



(22) 



is the dimensionless response function defined in the up- 
per half of the complex plane. The critical velocity is 



r 



defined as the velocity beyond which the amplitude of the 
collective excitations grows exponentially. In other words, 
if the imaginary part of any solution c changes its sign 
from negative to positive, then a transition from metasta- 
bility to instability occurs. The point at which this tran- 
sition occurs corresponds to the situation where Imc = 
in Eq. (|21|l and, consequently, the imaginary part in the 
r.h.s. of l|21(l must be equal to zero. This condition is 
satisfied provided that the imaginary part of Ea. l|22|l is 
zero which is the case only for Rec — 0. Thus we are left 
with a closed equation for the critical velocity v s that can 
be solved exactly: 



l k B T\ 3 th + 2gXn(0) 
k B TXf h - 2g Xn (0) 



lk B T\ 3 th -2gg 1/2 (e^-9no)) 
k B T\l h + 2gg 1/2 {e^-^)) 



(23) 
(24) 



It is common in the literature ( e.g. Q) to determine the 
critical velocity from the necessary condition that the re- 
lease energy of an excitation from the moving condensate 
becomes positive. In contrast, we use the stability condi- 
tion on the dynamics of the collective excitation which is 
necessary and sufficient to determine the critical velocity 
and this requires the knowledge of the Landau damping. 
As a consequence, this different criterion is not based on 
the energetic comparison between two thermodynamic 
states alone but rather on dynamical notions. 

The critical velocity is calculated as a function of the 
temperature using (|l(j[) . The dependences are displayed 
in Fig. 3 and 4 for gn/k B T c = 0.1 and 0.3 respectively. In 
both figures, the critical velocity is lower than the sound 
velocity for temperatures not too close to the transition 
point by an amount which can exceed 10 percents. This 
difference is due to the particular dependence of the self- 
energy matrix on the momentum and the frequency in 
the presence of depletion. This dependence is needed in 
order to calculate the sound velocity whereas only the 
static expressions are required for the critical velocity. 
For comparison, we have also plotted the sound veloc- 
ity obtained in the Popov approximation cp — y/ gn c /m 
1 1 1| . That approximation is gapless but not conserving 
and the self-energy matrix does not depend on k and u>. 
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Knowing that En = 2gn and £12 = gn c , the applica- 
tion of the previous reasoning to this approximation al- 
lows to deduce that the spectrum of collective excitations 
does not contain any imaginary part. Since no dynami- 
cal instability appears, only the condition of sign change 
in the energy spectrum allows to determine the critical 
velocity which in this case corresponds to the sound ve- 
locity. For T/T c > 0.8, we notice about a factor 2 of 
difference between the critical values obtained in the two 
approaches. As a consequence, the requirement that an 
approximation should be conserving modifies drastically 
the prediction. 
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FIG. 3: Critical velocity (full line) in units of y/gn/m as a 
function of the temperature for gn/ksT c = 0.1. The sound 
velocity (dashed line) and the Popov approximation cb (dot- 
ted line) are represented for comparison. 
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FIG. 4: Same as for FigOJbut for gn/k B T c = 0.3 



The influence of the relative velocity on the conden- 
sate population is also analyzed. Since the chemical po- 
tential depends on the relative velocity, the equation of 
state (|16|) shows that the condensate fraction decreases 
when the velocity increases. This result is plotted in Fig. 5 
where we compare the two situations without relative ve- 
locity and at critical velocity which show only a minor 
difference of about few percents. 

In conclusion, we determine the critical velocity from 
the dynamics of the collective excitations using the gap- 
less and conserving generalized random phase approx- 
imation. The calculated value displays relevant differ- 
ences in comparison with the Popov approach and ap- 
pears to be subsonic for a wide range of temperature 
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FIG. 5: Condensate fraction as a function of temperature for 
gn/kBTc = 0.3 at the critical velocity (full line) and in the 
absence of a relative velocity (dashed line). 



indicating that the presence of thermal excitations facil- 
itates the creation of an instability with the condensate. 
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